Transverse Entanglement Migration in Hilbert Space 
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We show that, ahhough the amount of mutual entanglement of photons propagating in free 
space is fixed, the type of correlations between the photons that determine the entanglement can 
dramatically change during propagation. We show that this amounts to a migration of entanglement 
in Hilbert space, rather than real space. For the case of spontaneous parametric down conversion, 
the migration of entanglement in transverse coordinates takes place from modulus to phase of the 
bi-photon state and back again. We propose an experiment to observe this migration in Hilbert 
space and to determine the full entanglement. 



Entanglement is one of the truly central features of 
the quantum world, and it forms the core of many ap- 
plications based on quantum theory. The observation 
of entanglement is generally achieved through the mea- 
surement of correlations between entangled subsystems. 
Correlation in quantum systems takes many forms and is 
open to observation in a variety of ways. Therefore, the 
determination of the amount of entanglement of quantum 
states depends on the measurement of the correlations 
where entanglement resides. This is of paramount im- 
portance, since in some experimental configurations one 
registers types of correlation that might not be appro- 
priate to quantify the entangled nature of the quantum 
state. 

In this Communication, we show that the measurement 
of correlation between paired photons can miss the detec- 
tion of entanglement entirely. The underlying reason is 
an interesting migration of entanglement that occurs in 
Hilbert space, but that depends on coordinate location 
in real space. This is manifest in photon correlations 
that show a rich and complex structure that evolves dur- 
ing propagation, although the amount of entanglement 
is constant. We focus here on entanglement that can be- 
come partly or entirely identified with the phase of the 
state, in which case the measurement of intensity correla- 
tions partially or completely misses the existing entangle- 
ment. This is an observable manifestation of the "phase 
entanglement" previously noted [l| for massive particle 
breakup in an Einstcin-Podolsky-Rosen (EPR) scenario. 

Entangled photons generated in spontaneous paramet- 
ric down-conversion (SPDC) are particularly open to the 
observation of this phenomenon. The generated two- 
photon states have been shown to exhibit entanglement 
in transverse momentum Q and in orbital angular mo- 
mentum [3, Moreover, one can enlarge the Hilbert 
space of the two-photon state by using several degrees of 
freedom 0] . The spatial transverse degrees of freedom of 
photon pairs produced in SPDC have attracted great at- 
tention because of the vast Hilbert space involved 0, 0] , 
and the availability of techniques to implement the d- 
dimensional quantum channel [1, [^, [l^l ■ 

Observations of SPDC entanglement have usually been 



made cither in the near zone or the far zone . Interest- 
ingly, in the course of photon propagation from the near 
field zone to the far field zone, the entanglement embed- 
ded in the two-photon positional amplitude migrates out 
of the positional wave function's modulus into its phase, 
and then back again. 

In the region between near and far zones, the entangle- 
ment not obtained through the measurement of intensity 
correlations can be recovered by measuring the phase in- 
formation of the joint wave function. Here we propose an 
experimental setup to accomplish this by exploiting the 
symmetries of the wave function. 

We consider a nonlinear optical crystal of length L, il- 
luminated by a quasi-monochromatic laser pump beam, 
propagating in the z direction. The signal and idler 
photons generated propagate from the output face of 
the nonlinear crystal under the sole effect of diffraction. 
The quantum state of the two-photon pair generated in 
SPDC, at a distance z from the output face of the non- 
linear crystal (z = 0), reads in wave number space as 
\^{z)) = J dp dq <^{p,q,z) a\{'p)a\{q)\Q,Q) , where p and 
q are the transverse wave numbers of the signal and idler 
photons, and aj(p) and a\{q) are the corresponding cre- 
ation operators. The signal and idler photons are as- 
sumed to be monochromatic. This assumption is justi- 
fied by the use of narrow band interference filters in front 
of the detectors. 

Under conditions of coUinear propagation of the pump, 
signal and idler photons with no Poynting vector walk- 
off, which would be the case of a noncritical type-II 
quasi-phase matched configuration, the mode function 
<I>(p, q,z) is given by 
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where iV is a normalization factor, A^. = kp{p + q) — 
ks{p) - h{cf} and Sk = kp{p + q) + ks{p) + h{q), Ep 
is the transverse profile of the pump at the input face 
of the nonlinear crystal, and kj {j = p,s,i) are wave 
number for the pump, signal, and idler waves. We have 
also made use of the paraxial approximation to describe 
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the propagation of the signal and idler photons in free 
space. 

The sine function that appears in Eq. ([1]) can be ap- 
proximated by a Gaussian exponential function which 
accurately retains the main features of the entanglement 
of the wave function We take sincte^ ~ cxp[— aba:^] 
with a = 0.455, so that both functions coincide at the 
1/e^ intensity. Here we assume a pump beam with a 
Gaussian shape. Therefore, the mode function can be 
written as [12, llJ 
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where N = {{wlaL)/[TT^l + w*/<jl)k^]Y^\ fi,{z) = 
2{z+L)/ep~ao/il+al/wt) and/X2(z) ^ {2z+L)/k°. We 
denote wq as the pump beam width and ctq = ^'^R/kp, 
with R being the radius of curvature of the Gaussian 
beam at the entrance face of the nonlinear crystal and 
kp = ajpHp/c. ujp and Up are the corresponding angular 
frequency and refractive index, respectively. Notice that 
we have made use of the approximation kp ~ 2fc" = 2fc". 
Moreover, all phase factors independent of the transverse 
variables have been neglected. 

Equation ^ describes the two-photon quantum 
state in transverse wave number space (p, g). We 
can also describe the two-photon quantum state in 
coordinate space. In this case, "^{xs^Xi,z) = 
l/(27r)^ j dp dq <E>(p, q, z) exp(ip ■ Xg + iq- Xi), and since 
Eq. ([2) can be written as ^{p,q,z) = F{p + q,z)G{p — 
q, z), one can easily obtain 
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where N' ^ N/{-iPY'^., (3{z) = {aL/klf + ^4{z), and 
7(2;) = ^0/(1 + ifo/o-g)'^ + The conditional coinci- 

dence rate in coordinate space is given by TZxi^s, Xi, z) = 
|4'(xs, while the conditional coincidence rate in 

momentum space is 7^p(p, q, z) = \^(p,q,z)\'^. 

Equation ^ shows that the two-photon state is sep- 
arable in momentum space, and in coordinate space, if 
the conditions Wq/{1 + Wq/o-q) ~ aL/kp and (To/(1 + 
(Tq/wq) = L/kp are fulfilled, which corresponds to sep- 
arability in modulus and phase. Notwithstanding, from 
Eq. ([3]) we also observe that it is possible that the bi- 
photon function is separable in modulus at a specific loca- 
tion, although not in phase. Therefore, |^'(afs, a?;, ^o)P = 

\'^siXs,Zo)\^\■^^ix^,z)\^. 

The central point is that, at a certain location zq from 
the output face of the nonlinear crystal, where the z- 
dependent condition kp'WQf3{zo) = {1 + Wq / (jQ)aL^{zo) is 
fulfilled, one does not observe intensity correlations at all 
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FIG. 1; (a) Ellipticity of the bi-photon function in coordinate 
space, (b) Fedorov ratio J-x in coordinate space. The dashed 
line corresponds to TZx ~ 1- (c) The spatial conditional co- 
incidence rate TZx{xs, Xi, z) at different locations: z — Ocm, 
z = Zq — 6.2 cm and at z — 20 cm. Parameters: Crystal 
length L — 5 mm; pump beam width wo ~ 800 /im; pump 
beam wavelength Ap = 800 nm. The ellipticity is plotted in 
logarithmic scale. 



in coordinate space, although the quantum state is not 
separable in either momentum or coordinate. Since the 
amount of entanglement is determined by the existing 
correlations of the bi-photon function in modulus and 
phase, at zq all entanglement lives in the phase of the 
bi-photon function in coordinate space. 

Figure[lja) shows the evolution of the intensity correla- 
tions as a function of the distance z. We plot the elliptic- 
ity (e) in the plane {xs+Xi, Xs — Xi) of the bi-photon func- 
tion given by Eq. ([3]), i.e., e = k^w^P/[{l +WQ/aQ)aLj]. 
The spatial conditional coincidence rate TZ^ {xs ,Xi,z) for 
three specific locations z are also shown in the figure. 
Note that the amount of entanglement does not depend 
on the location z, while the magnitude of the intensity 
correlations evolves with z, as shown by the variation 
of ellipticity in Fig. [TJ Therefore, although the amount 
of entanglement is unchanged with z, the type of corre- 
lation that determines the entanglement is different at 
every location z. 

In order to quantify the amount of entanglement of 
the two-photon state, we perform the Schmidt decompo- 
sition [ij, EBl of the bi-photon function given by Eq. ^ 
at the output face of the nonlinear crystal. As shown 
in Appendix A, the amount of entanglement denoted K 
(i.e., the Schmidt number for continua is given by 



K 



[^{A + B)f + [Q{A-B)f 



(4) 



where A = Wq/{1 + Wq/gq) + i^i and B = aL/kp + 1^2- 
Note that K does not depend on z even though /ii and 
fi2 do so. 

For the sake of comparison, let us consider the Fe- 
dorov ratio [1^, here denoted JF, a typical correla- 
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FIG. 2: (Color online) Experimental scheme to detect the to- 
tal entanglement. The signal photon is directed into a modi- 
fied Mach-Zehnder interferometer with two Dove prisms DPi 
with orientation angles 9i. PBS is a polarization beam split- 
ter. 

tion measurement that could potentially be employed 
to show the existence of entanglement. For the signal 
photon in momentum space it takes the form Ts^p = 
(A'^Ps) / {A'^'p's)i, and the expression in coordinate space 
is Ts^x = {A'^Xs) / {A'^Xs)i. Here the variance averages 
not containing subscript i arc unconditional. The aver- 
ages with subscript i are conditioned on the idler photon, 
which is to be constrained by pi = and = 0, in the 
Ps and Xs averages respectively. 

If the entanglement resides only in the modulus of the 
bi-photon given by Eq. i.e., fii — fj,2, can be 
shown to be equal to the amount of entanglement given 
by Eq. (j4]), while J-x only gives the correct amount of 
entanglement in the near and far fields. This is the typ- 
ical experimental condition if the pump beam shows no 
curvature at the input face of the nonlinear crystal. 

If part or all the entanglement resides in the phase of 
the bi-photon. even Tp does not correctly measure the 
amount of the entanglement of the quantum state, only 
the part of the entanglement that resides in the modulus 
of the bi-photon function. Figure[T]^b) shows the Fedorov 
ratio in coordinate space for a typical case. At z = zqj 
where the bi-photon function shows no cUipticity in the 
modulus, we have Tx = 1, although the quantum state 
is entangled. 

In Fig.[2]wc show an experimental scheme to detect the 
total entanglement of the bi-photon described by Eqs. ([2|) 
or The signal photon is sent to a modified Mach- 
Zchnder interferometer with two Dove prisms inserted 
in the interfering arms. The arms are assumed to be 
balanced so that the relative phase shift between the two 
arms of the interferometer due to propagation is zero. 

We set the orientation angles of the Dove prisms 9i ^ 
tt/2 and 62 = 0. The conditional coincidence rates of the 
output ports of the interferometer shown in Fig. [2] take 
the form (see Appendix B) 
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FIG. 3: (Color online) The coincidence interference pattern 
observed in the setup of Fig. (2] The dotted line shows the 
situation when phase entanglement is absent. 

where 

1 2 

Pa,,ai{Xs,Xi) = - \'if{Xs,ys]X,) + * (-Xs , - J/s ; .f ) | , (6a) 
1 2 

Pb,,a,{xs,Xi) = - \'if{xs,-ys;xi) - ^'(-Xs, y^; ^j)| ■ (6b) 

Therefore, the amount of entanglement of the quantum 
state given by Eq. ^ can be quantified as K = (P+ + 
P_ ) / (P-f — P_ ) . The experimental setup plotted in Fig. [2] 
measures the "full" entanglement of the quantum state 
described by Eqs. ©or p. 

On the other hand, the joint probability distributions 
in Eq. ([5]) also exhibit interesting interference behavior. 
Using Eq. ([3]), we have 

Pdifii^Xg ; X^ ) — Pas .ai {Xg , X^) Pbs ,ai (^s ; Xi) 

= AA2e-^+(")(^'+^-) cos [21^{z)xs ■ x,] , (7) 

where i?+(z) ^ wH [{I + w^/ al)^{z)\ + aL/ [fcO/3(z)] 
and /-(z) = /ii(z)/7(z) — /i2(z)//3(z). The diff'erence of 
the joint probability Pdiff is plotted in Fig.[3]as a function 
of Xs ■ It is seen that the location of the second maximum 
is at Xs ~ 0.9537r/a;i/_, from which we can determine the 
entanglement in the phase. That is, the interference of 
the wave function with itself, with its symmetry, reveals 
not only the phase information of the wave function, but 
more importantly how the two photons are correlated in 
the phase through the fringe pattern given by Eq. ([7|). 

In conclusion, we have demonstrated the Hilbert space 
migration of entanglement of down-converted photons in 
free-space propagation. We suggested its implication for 
experiments involving the quantification of the degree of 
entanglement by means of common variance measure- 
ments. We have also suggested a simple experimental 
scheme that can detect both the entanglement in ampli- 
tude and phase. 

KWC acknowledges support from a fellowship from 
the Croucher Foundation; JPT from the Generalitat de 
Catalunya, from the European Commission under the 
integrated project Qubit Applications (QAP) funded by 
the 1ST directorate (Contract No. 015848), and from the 
MEC (Consohdcr Ingcnio 2010 project 0/OrCSD2006- 
00019, FIS2004-03556); and JHE from ARO Grant 
W911NF-05-1-0543. Discussions with R.W. Boyd, C.K. 
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Law, M.N. O' Sullivan-Hale and K. Wodkicwicz have been 
very helpful. 

Appendix A — Let us consider a bi-photon amplitude 
that, at the output face of the nonlinear crystal (z = 0), 
is written as 



3?(A)5R(B)ii 



exp 



A\p+\' 



B\p-\ 



, (8) 



where = P + Q and p_ = p — q. Inspec- 

tion of Eq. ([5]) shows that the bi-photon function can 
be separated for the two transverse dimensions, i.e., 

y,qy)- Therefore, the Schmidt 
decomposition of Eq. ([5]) can be written as g) = 

Sm.n^O V Amn/mn {p)g7nn (<z) ; whcrC the basis fuUCtiouS 

of the decomposition are fmnip) = 4'm{px)'4'n{Py) and 
gmn{q) = '>Pmiqx)tpniqy) , with eigenvalue VAmA„. 

The reduced density matrix for the signal pho- 
ton, ps(j>x,Py,Px,Py,) = Tridior I Can be sep- 
arated into two matrices , i.e., ps {px ,Py,Px,Py,) = 
Pxipx , Px) PyiPy , Py , ), each One corresponding to a trans- 
verse coordinate. The functions and the eigen- 
values Xm can be obtained from the expansion of the 
one dimensional reduced density matrix px in the form 

Px{Px,Px) = Y.'^=0'^^.ri'ipn{Px)^n{Px)- 

From Eq. ^ , the one-dimensional reduced density ma- 
trix becomes 



PxiPx,Px) 



l2^{A)^{B) 



■ exp 



y 7r3?(A + B) ^ t 2^{A + B) 
X exp {-[(a + b){pl+pl) - 2bpxPx]} , (9) 



where a = n{A)n{B) /di{A + B) and b ^ \A - 
B\'^/[8'Si{A + B)]. The representation given by Eq. ^ 
is also found in the determination of the mode structure 
of Gaussian-Schell model sources in the theory of partial 
coherence One can find that the eigenvalue Xx.n is 
given by Xx,n ~ (a/c)^/^(l — ti;^)^/^it;", with c = (a^ + 
2a6)^/^ and w = b/ (a+&+c). The functions tpn are appro- 
priately normalized Hermite-Gaussian functions. There- 
after, the calculation K = 1/ X^m n=o('^2:,mAy,„)^ 



would yield K = (c/a)^, which can be straightforwardly 
written as the expression that appears in Eq. @ . Notice 
that K describes entanglement in the {px,Py) space. 

Appendix B — The modified Mach-Zehnder interfer- 
ometer shown in Fig. [2] contains three basic elements 
that modify the spatial shape of the bi-photon function. 
The action of the mirrors is described by ain{x,y) 
o.out{x, —y), where x and y are the transverse coordinates 
in the frame of each individual beam. 

The action of the beam-splitter is ain{x,y) 
at{x,y) + iar{x, —y), where at is the creation operator 
of the transmitted photon, and the corresponding 
creation operator of the reflected photon. For a Dove 
prism that is rotated by an angle 6 with respect to 
the axis of image inversion, the fields before and after 
the dove prism are given by ain[x,y) — > aout{xcos29 — 
y sin 29,— X sin 26 — y cos 29) . Together with the effect of 
the polarization beam-splitter, one thus obtains that all 
joint probability detections in the configuration described 
in Fig. [5] are given by Eq. 

The bi-photon function at location z can be written in 
the form 

M 1/2 



a{l-w^) 



E 

n.m—O 



x-0„i(xs, z)^n{ys,z)ipm{xi, z)il;„{yi, z) , (10) 

where the function ipn{x,z) corresponds to Hermite- 
Gaussian function at z = that evolves under the 
sole influence of diffraction. Due to the symmetry of 
the Hermite-Gaussian functions, one has x, z) = 

(— l)"'!/)„(a;, z). Making use of this symmetry property, 
one obtains 



P+ = 



a{l-w^) 
2^ 



n.m—O 



7i-\-'m 



+ i-w) 



(11) 



From Eq. ((TT|), one obtains Eq. ([5]), taking into account 
that the amount of entanglement is given by K = (cj a)"^ . 
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